Optomechanical detection of weak forces 



David Vitali, Stefano Mancini and Paolo Tombesi 
INFM and Dept. of Physics, University of Camerino, 1-62032 Camerino Italy 

ABSTRACT 

Optomechanical systems are often used for the measurement of weak forces. Feedback loops can be used in these 
systems for achieving noise reduction. Here we show that even though feedback is not able to improve the signal 
to noise ratio of the device in stationary conditions, it is possible to design a nonstationary strategy able to 
' improve the sensitivity. 
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1. INTRODUCTION 

' Optomechanical devices are used in high sensitivity measurements, as the interferometric detection of gravita- 
tional waves, 1 and in atomic force microscopes. 2 Up to now, the major limitation to the implementation of 
t-H sensitive optical measurements is given by thermal noise. 3 Some years ago it has been proposed 4 to reduce 
thermal noise by means of a feedback loop based on homodyning the light reflected by the oscillator, playing 
the role of a cavity mirror. This proposal has been then experimentally realized 5-7 using the "cold damping" 
technique, 8 which is physically analogous to that proposed in Ref. 4 and which amounts to applying a viscous 
feedback force to the oscillating mirror. In these experiments, the viscous force is provided by the radiation 
pressure of another laser beam, intensity-modulated by the time derivative of the homodyne signal. 

Both the scheme of Ref. 4 and the cold damping scheme of Refs. 5 ~ 7 cool the mirror by overdamping it, 
thereby strongly decreasing its mechanical susceptibility at resonance. As a consequence, the oscillator does not 
resonantly respond to the thermal noise, yielding in this way an almost complete suppression of the resonance peak 
in the noise power spectrum, which is equivalent to cooling. However, the two feedback schemes cannot be directly 
applied to improve the detection of weak forces. In fact the strong reduction of the mechanical susceptibility at 
resonance means that the mirror does not respond not only to the noise but also to the signal. This means that 
the signal to noise ratio (SNR) of the device in stationary conditions is actually never improved. 9 ' 10 Despite 
that, here we show how it is possible to design a nonstationary strategy able to significantly increase the SNR for 
the detection of impulsive classical forces acting on the oscillator. This may be useful for microelectromechanical 
systems, where the search for quantum effects in mechanical systems is very active, 11 as well as for the detection 
of gravitational waves. 1 
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2. THE MODEL 

We shall consider a simple example of optomechanical system, a Fabry-Perot cavity with a movable end mirror 
(see Fig. 1 for a schematic description) . The optomechanical coupling between the mirror and the cavity field is 
realized by the radiation pressure. The electromagnetic field exerts a force which is proportional to the intensity 
of the field, which, at the same time, is phase-shifted by an amount proportional to the mirror displacement from 
its equilibrium position. The mirror motion is the result of the excitation of many vibrational modes, including 
internal acoustic modes. However here we shall study the spectral detection of weak forces acting on the mirror 
and in this case one can focus on a single mechanical mode only, by considering a detection bandwidth including 
a single mechanical resonance peak. 12 Optomechanical devices are able to reach the sensitivity limits imposed 
by quantum mechanics 10 and therefore we describe the mechanical mode as a single quantum harmonic oscillator 
with mass m and frequency oj m . 

When iv m is much smaller than the cavity free spectral range, one can focus on one cavity mode only, because 
photon scattering into other modes can be neglected. 13 Moreover, in this regime, the generation of photons due 




Figure 1. Schematic description of a linear Fabry-Perot cavity with the end oscillating mirror M. The equilibrium cavity 
length is L. A cavity mode is driven by an input laser beam. The output field is subjected to homodyne detection (D). 
The signal is then fed back to the mirror motion (dashed line). 

to the Casimir effect, and also retardation and Doppler effects are completely negligible. 14 The system is then 
described by the following two-mode Hamiltonian 14 

H = hujctfb + huj m (P 2 + Q 2 ) - 2hGtfbQ - 2hf(t)Q + ihE (&t e -*"°* - be iuJot ) , (1) 

where b is the cavity mode annihilation operator with optical frequency uj Ci f(t) is the classical force to be 
detected, and E describes the coherent input field with frequency ujq ~ to c driving the cavity. The quantity 
E is related to the input laser power p by E = y pj c /ft,u>o, being j c the photon decay rate. Moreover, Q 
and P are the dimensionless position and momentum operator of the movable mirror. It is [Q,P] = i/2, and 
G = (lo c / L)^Jh/2muj m represents the optomechanical coupling constant, with L the equilibrium cavity length. 

The dynamics of an optomechanical system is also influenced by the dissipative interaction with external 
degrees of freedom. The cavity mode is damped due to the photon leakage through the mirrors which couple 
the cavity mode with the continuum of the outside electromagnetic modes. For simplicity we assume that 
the movable mirror has perfect reflectivity and that transmission takes place through the fixed mirror only. 
The mechanical mode undergoes Brownian motion caused by the uncontrolled coupling with other internal and 
external modes at thermal equilibrium. We shall neglect in our treatment all the technical noise sources: we 
assume that the driving laser is stabilized in intensity and frequency and we neglect the electronic noise in 
the detection circuit. Including these supplementary noise sources is however quite straightforward (see for 
example 15 ). Moreover recent experiments have shown that classical laser noise can be made negligible in the 
relevant frequency range. 3, 16 

The dynamics of the system can be described by the following set of coupled quantum Langevin equations 
(QLE) (in the interaction picture with respect to tiuiowb) 

Q(t) = uj m P(t), (2) 
P(t) = -Lu m Q(t)+W(t) + f(t)~ lm P(t)+Gb\t)b(t) 7 (3) 

b(t) = -(iLu c -iu + ^b(t) + 2iGQ{t)b(t)+E + ^T c b m (t), (4) 

where bi n (t) is the input noise operator 1 associated with the vacuum fluctuations of the continuum of modes 
outside the cavity, having the following correlation functions 

(b in (t)b in {t')) = (b\ n (t)b m (t')) = , (5) 
(b in (t)bl(t')) = 6(t-t'). (6) 

Furthermore, W(i) is the quantum Langevin force acting on the mirror, with the following correlation function, 18 



(7) 



with T the bath temperature, 7 m the mechanical decay rate, ks the Boltzmann constant, and w the frequency 
cutoff of the reservoir spectrum. 

In standard applications the driving field is very intense so that the system is characterized by a semiclassical 
steady state with the internal cavity mode in a coherent state \(3), and a new equilibrium position for the mirror, 
displaced by G\(i\ 2 /u m with respect to that with no driving field. The steady state amplitude (3 is given by the 
solution of the classical nonlinear equation (3 = E/ (j c /2 + ioj c — iojo + 2iG 2 /uj m \{3\ 2 ). In this case, the dynamics is 
well described by linearizing the QLE (2)-(4) around the steady state. Redefining with Q(t) and b(t) the quantum 
fluctuations around the classical steady state, introducing the field phase Y(t) = i (b^(t) — b(t)) /2 and field 
amplitude X(t) = (&(i) + b^(t)) /2, and choosing the resulting cavity mode detuning A = uj c — uj +2G 2 / Lo m j3 2 = 
(by properly tuning the driving field frequency w )j the linearized QLEs can be rewritten as 

Q(t) = w m P(t), (8) 
P(t) = -uj m Q(t)- lm P(t) + 2G[3X(t)+W(t)+f(t), (9) 

Y(t) = -^Y(t) + 2G/3Q(t) + ^Y in (t), (10) 
X{t) = -^X(t) + ^X in (t), (11) 

where we have introduced the phase input noise Y in (t) — i (b\ n (t) — b in (t)j and the amplitude input noise 
X in {t) = b\ n {t) + b in {t). 

3. POSITION MEASUREMENT AND FEEDBACK 

As it is shown by Eq. (10), when the driving and the cavity fields are resonant, the dynamics is simpler be- 
cause only the phase quadrature Y(t) is affected by the mirror position fluctuations Q(t), while the amplitude 
quadrature X(t) is not. Therefore the mechanical motion of the mirror can be detected by monitoring the phase 
quadrature Y{t). The mirror position measurement is commonly performed in the large cavity bandwidth limit 
7 C 3> Gj3, Ld m , when the cavity mode dynamics adiabatically follows that of the movable mirror and it can be 
eliminated, that is, from Eq. (10), 

Y(t)c^l Q (t) + ^B, (12) 

and X(t) ~ X in (t) / \/^ from Eq. (11). The experimentally detected quantity is the output homodync photocur- 
rent 19 - 21 

Y out {t) = 2 v ^Y(t) - y/rjY? n (t) , (13) 

where -q is the detection efficiency and Y^ n {t) is a generalized phase input noise, coinciding with the input noise 
Yin{t) in the case of perfect detection 77 = 1, and taking into account the additional noise due to the inefficient 
detection in the general case r\ < l. 20 This generalized phase input noise can be written in terms of a generalized 
input noise b v (t) as Y^ n {t) = i [b^(t) — b v (t)\ . The quantum noise b v (t) is correlated with the input noise 6j n (i) 
and it is characterized by the following correlation functions 20 

(b r ,(t)b ri (t')) = (bl 1 (t)b ri (t'))=0, (14) 
(b ri (t)bl l (t'))=S(t-t'), (15) 
<& <n (t)&t (f )> - (br,(t)b\ n (t')) = VvS(t f )■ (16) 

The output of the homodyne measurement may be used to devise a phase-sensitive feedback loop to control the 
dynamics of the mirror, as in the original proposal, 4 or in cold damping schemes. 5-8 Let us now see how these 
two feedback schemes modify the quantum dynamics of the mirror. 

In the scheme of Ref., 4 the feedback loop induces a continuous position shift controlled by the output 
homodyne photocurrent Y out (t). This effect of feedback manifests itself in an additional term in the QLE for a 



generic operator 0(t) given by 



Of b {t) = ( dt'G mf (t')Y out (t - t 1 ) [g mf P(t), 0{t)] , (17) 
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where G m f(t) is the feedback transfer function, and g m f is a feedback gain factor. The implementation of this 
scheme is nontrivial because it is equivalent to add a feedback interaction linear in the mirror momentum, as it 
could be obtained with a charged mirror in a homogeneous magnetic field. For this reason here we shall refer to 
it as "momentum feedback" (see, however, the recent parametric cooling scheme demonstrated in Ref., 7 showing 
some similarity with the feedback scheme of Ref. 4 ). 

Feedback is characterized by a delay time which is essentially determined by the electronics and is always 
much smaller than the typical timescale of the mirror dynamics. It is therefore common to consider the zero 
delay-time limit G m f(t) ~ S(t). For linearized systems, the limit can be taken directly in Eq. (17), 10 so to get 
the following QLE in the presence of feedback 

Qit) = uj m P(t)+g mnc Y(t)- g -^ J^Y? n (t), (18) 

z V " 

P(t) = -cj m Q(t)- 7m P(t) + 2GpX(t)+W(t) + f(t), (19) 

Y(t) = -^Y(t)+2Gf3Q(t) + ^Y in (t), (20) 

X(t) = -^X(t) + ^X in (t), (21) 

where we have used Eq. (13). After the adiabatic elimination of the radiation mode (see Eq. (12)), and introducing 
the rescaled, dimensionless, input power of the driving laser £ = 16G 2 /3 2 /7 m 7 c = 64G 2 p/?iu;o7m7c , and the 
rescaled feedback gain gi = —AGl3g m f /7 TO , the above equations reduce to 

Q(t) = WmP (t)- lm g 1 Q{t)- s p^g 1 Y in {t) + ^ 9 -^Yl> n {t), (22) 



P(t) = -io m Q(t)- lm P(t) + -^CX m (t)+W(t) + f(t). (23) 

This treatment explicitly includes the limitations due to the quantum efficiency of the detection, but neglects 
other possible technical imperfections of the feedback loop, as for example the electronic noise of the feedback 
loop, whose effects have been discussed in. 6 

Cold damping techniques have been applied in classical electromechanical systems for many years, 8 and only 
recently they have been proposed to improve cooling and sensitivity at the quantum level. 22 This technique is 
based on the application of a negative derivative feedback, which increases the damping of the system without 
correspondingly increasing the thermal noise. 22 This technique has been succesfully applied to an optomechanical 
system composed of a high- finesse cavity with a movable mirror in. 5 ~ 7 In these experiments, the displacement 
of the mirror is measured with very high sensitivity, 3 ' 7 and the obtained information is fed back to the mirror 
via the radiation pressure of another, intensity-modulated, laser beam, incident on the back of the mirror. Cold 
damping is obtained by modulating with the time derivative of the homodyne signal, in such a way that the 
radiation pressure force is proportional to the mirror velocity. A quantum description of cold damping can 
be obtained using either quantum network theory, 22 or a quantum Langevin description. 9, 10 In this latter 
treatment, cold damping implies the following additional term in the QLE for a generic operator 0(t), 

6 fb (t) = ^— f dt'G cd (t')Y out (t - t') [g cd Q(t),0(t)] , (24) 
V\TU Jo 



where G c d{t) and g cd are the corresponding transfer function and gain factor. As in the previous case, one 
usually assume a Markovian feedback loop with negligible delay. Since one needs a derivative feedback, this 



would ideally imply G c d(t) = —5'(t), i.e., G c d{oj) = iu>, Vw, even though, in practice, it is sufficient to satisfy this 
condition within the detection bandwidth. In this case, the QLEs for the cold damping feedback scheme become 



Q(t) = w m P(t), (25) 

led -, 



P(t) = -u> m Q(t) - j m P(t) + 2G(3X{t) - g cd Y(t) + -M^Y?Jt) + W{t) + f(t) , (26) 



Y(t) = -^ Y (t) + 2G(3Q(t) + ^Y m (t), (27) 



X(t) = -^ X (t) + ^X in (t). (28) 

Adiabatically eliminating the cavity mode, and introducing the rescaled, dimensionless feedback gain g 2 = 
4G/3w m # cd /7 m 7 c , one has 

Q(t)=u m P(t), 

P(t) = -u m Q(t) 7ra P(t) + ^X m (t) + W(t) + f{t) - ^Q(t) - 9 -^Y m (t) + f^YHt). 

The presence of an ideal derivative feedback implies the introduction of two new quantum input noises, Y in (t) 
and Y^l(t), whose correlation functions can be simply obtained by differentiating the corresponding correlation 
functions of Y in (t) and Y^(t). 10 However, as discussed above, these "differentiated" correlation functions have 
to be considered as approximate expressions valid within the detection bandwidth only. 

The two sets of QLE for the mirror Heisenberg operators show that the two feedback schemes are not exactly 
equivalent. They are however physically analogous, as it can be seen, for example, by looking at the Fourier 
transforms of the corresponding mechanical susceptibilities 9, 10 

Xm/M = uj m [u4 + 017^ - u? + iuj lm (1 + gi )] _1 (29) 
for the momentum feedback scheme, and 

XcdM = Wm [uf n - LO 2 + iw-fm (1 +52)] 1 (30) 

for cold damping. These expressions show that in both schemes the main effect of feedback is the modification of 
mechanical damping -f m — > 7m(l + 9i) {i — 1, 2). Therefore, also momentum feedback provides a cold damping 
effect of increased damping without an increased temperature. In this latter case, one has also a frequency 
renormalization lu^ — > + 7„<?i, which is however negligible when the mechanical quality factor Q = u m /jm 
is large. 



4. SPECTRAL MEASUREMENTS AND THEIR SENSITIVITY 



Spectral measurements are performed whenever the classical force f(t) to detect has a characteristic frequency. 
We adopt a very general treatment which can be applied even in the case of nonstationary measurements. The 
explicitly measured quantity is the output homodyne photocurrent Y out (t), and therefore we define the signal 
S(uj) as 
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dte- wt (Y out (t))F Tm (t) 



(31) 



where Fr m (t) is a "filter" function, approximately equal to one in the time interval [0, T m ] in which the measure- 
ment is performed, and equal to zero otherwise. The spectral measurement is stationary when T m is very large, 
i.e. is much larger than all the typical timescales of the system. Using Eq. (12) and the input-output relation 
(13), the signal can be rewritten as 



S(w) 
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ckj'x(u>')f(w')F Tm (u>-u>') 



(32) 



where f(u>) and Fr m {uj) are the Fourier transforms of the force and of the filter function, respectively, and x(u>) 
is equal to x m /(w) or Xcd(w), according to the feedback scheme considered. 

The noise corresponding to the signal S(ui) is given by its "variance" ; since the signal is zero when f(t) = 0, 
the noise spectrum can be generally written as 



N(u) 



a +00 r+00 
dtF Tm (t) / dt'FT ra (t')e- <w (*- t ')<y out (t)y out (t')) /= o 
00 J —00 



1/2 



(33) 



where the subscript / = means evaluation in the absence of the external force. Using again (12), Eqs. (13), 
and the input noises correlation functions (5)-(6) and (14)-(16), the spectral noise can be rewritten as 



JVM = ( (8G ^ )2 / +OC dtF Tm (t) f + °° dt'F Tm (t')e 
I 7c J — 00 J — 00 



-iuj{t — t' 



] C{t,t') 



+ 00 -\ !/ 2 

dfcFr m (t) 2 l , (34) 



where C(t,t') — (Q(t)Q(t') + Q(t')Q(t)) /2 is the symmetrized correlation function of the oscillator position. 
This very general expression of the noise spectrum is nonstationary because it depends upon the nonstationary 
correlation function C(t, t'). The last term in Eq. (34) is the shot noise term due to the radiation input noise. 

4.1. Stationary spectral measurements 

Spectral measurements are usually performed in the stationary case, that is, using a measurement time T m 
much larger than the typical oscillator timescales. The most significant timescale is the mechanical relaxation 
time, which is 7" 1 in the absence of feedback and [7 m (l + ffi)] -1 {i = 1,2) in the presence of feedback. In 
the stationary case, the oscillator is relaxed to equilibrium and, redefining t' = t + r, the correlation function 
C(t, t') = C(t,t + r) in Eq. (34) is replaced by the stationary correlation function C s t{r) — lim^oo C(t, t + r). 
Moreover, for very large T m , one has Fr m (t + r) ~ FT m {t) ~ 1 and, defining the measurement time T m so that 
T m = J dtF Tm (t) 2 , Eq. (34) assumes the form 



N(oj) 



■N 2 Q {uj) + n 
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where 



+00 



dre-^ T C(r), 



(35) 



(36) 



is the stationary position noise spectrum. This noise spectrum can be evaluated by solving the quantum Langevin 
equations in the presence of the two feedback schemes and Fourier transforming . 10 One obtains 
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for the momentum feedback scheme and 



2u n 
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(37) 



(38) 



for the cold damping scheme, where Gi(u) (i=mf,cd) are the Fourier transforms of the feedback transfer functions 
and &i(u)) is a gate function equal to 1 within the frequency interval / and zero outside. The position noise 
spectrum for the momentum feedback essentially coincides with that already obtained in, 4 except that in that 
paper the high-temperature (coth^w/^fceT) ~ 2ksT /hu), Markovian feedback (G m /(u>) ~ 1), and infinite cutoff 
(m — > 00) approximations have been considered. The noise spectrum in the cold damping case of Eq. (38) instead 
essentially reproduces the one obtained in, 23 with the difference that in Ref. 23 the homodyne detection efficiency 
r\ is set equal to one, and the Markovian feedback (G c d(u>) — iui), and infinite cutoff (w — > 00) approximations 
have been again considered. The comparison between Eqs. (38) and (37) shows once again the similarities of 
the two schemes. The only differences lie in the different susceptibilities and in the feedback-induced noise term, 



which has an additional factor in the momentum feedback case, which is however usually negligible with 

good mechanical quality factors. In fact, the two noise spectra are practically indistinguishable in a very large 
parameter region. 

The effectively detected position noise spectrum is not given by Eqs. (38) and (37), but it is the noise 
spectrum associated to the output homodyne photocurrent of Eq. (35) rescaled to a position spectrum. This 
homodyne-detected position noise spectrum is actually subject also to cavity filtering, yielding an experimental 
high frequency cutoff 7 C , which however does not appear in our expressions because we have adiabatically 
eliminated the cavity mode from the beginning. Therefore the noise spectrum derived here is correct only for 
u> < 7 C . However, this is not a problem because the frequencies of interest (i.e. those within the detection 
bandwidth) are always much smaller than j c , and also than w. Within the detection bandwidth one can safely 
approximate G m f(u) ~ I, G c d(oj) — iu and 6[_ TO;IJ7 ](u;) = 1 in Eqs. (38) and (37). Therefore, using Eq. (35), 
the detected position noise spectrum can be written 



N Q,det(u) = 7m fcMI 2 



+ coth 



4??C7 



1 (39) 
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where i — 1 refers to the momentum feedback case and i = 2 to the cold damping case. 

This spectrum has four contributions: the radiation pressure noise term, proportional to the input power 
p, the feedback-induced term proportional to the squared gain and inversely proportional to p, the Brownian 
motion term which is independent of p, and the shot noise term inversely proportional to p. The main effect of 
feedback on the noise spectrum is the modification of the susceptibility due to the increase of damping, yielding 
the suppression and widening of the resonance peak. This peak suppression in the noise spectrum has been 
already predicted and illustrated in, 4, 23 and experimentally verified for the cold damping case in. 5, 6 It has 
been shown 10 ' 23 that both feedback schemes are able to arbitrarily reduce the displacement noise at resonance. 
This noise reduction at resonance is similar to that occurring to an oscillator with increasing damping, except 
that in the present case, also the feedback-induced noise increases with the gain, and it can be kept small only 
if the input power is correspondingly increased in order to maintain the optimal value minimizing the noise. 10 
This arbitrary reduction of the position noise in a given frequency bandwidth with increasing feedback gain does 
not hold if the input power £ is kept fixed. In this latter case, the noise has a frequency-dependent lower bound 
which cannot be overcome by increasing the gain. 

In the case of stationary spectral measurements also the expression of the signal simplifies. In fact, one has 
F Tm {uj) ~ 5{u), and Eq. (32) becomes S(u>) = 8G(3r]\x{uj)f(u>)\/2Tv^%. The stationary SNR, 1Z st (w), is now 
simply obtained dividing this signal by the noise of Eq. (35), 



n st (co) = |/>)K 7m T„ 



2u m \2k B T) 4 4r,CU^ v ^ > m > ^ m \U"f 



(40) 



where again i — 1 refers to the momentum feedback case and i = 2 to the cold damping case. It is easy to see 
that, in both cases, feedback always lowers the stationary SNR at any frequency, (except at w = 0, where the 
SNR for the cold damping case does not depend upon the feedback gain). This is shown in Fig. 2, where the 
stationary SNR in the case of an ideal impulsive force (that is, f(u>) is a constant) is plotted for three values 
of the feedback gain. The curves refer to both feedback schemes because the two cases i = 1, 2 gives always 
practically indistinguishable results, except for very low values of Q. This result can be easily explained. In 
fact, the main effect of feedback is to decrease the mechanical susceptibility at resonance (see Eqs. (29) and 
(30)), so that the oscillator is less sensitive not only to the noise but also to the signal. Therefore, even though 
the two feedback schemes are able to provide efficient cooling and noise reduction in narrow bandwidths for the 
mechanical mode, they cannot be used to improve the sensitivity of the optomechanical device for stationary 
measurements. In the next section we shall see how cooling via feedback can be used to improve the sensitivity 
for the detection of impulsive forces, using an appropriate nonstationary strategy. 

5. HIGH-SENSITIVE NONSTATIONARY MEASUREMENTS 



The two feedback schemes discussed here achieve noise reduction through a modification of the mechanical 
susceptibility. However, this modification does not translate into a sensitivity improvement because at the same 




Figure 2. Stationary SNR as a function of frequency in the case of an ideal impulsive force, i.e., / (o>) = const. The 
full line refers to the case with no feedback, the dashed line to the case with g 1 = g 2 = 10 4 , and the dotted line to the 
case with g\ = gi = 10 5 (the two feedback schemes give indistinguishable results in these cases). The other parameters 
are Q = 10 5 , £ = 10, ksT '/ftiOm = 10 5 , and r\ = 0.8. At a given frequency, the stationary SNR decreases for increasing 
feedback gain. 

time it strongly degrades the detection of the signal. The sensitivity of position measurements would be improved 
if the oscillator mode could keep its intrinsic susceptibility, unmodified by feedback, together with the reduced 
noise achieved by the feedback loop. This is obviously impossible in stationary conditions, but a situation very 
similar to this ideal one can be realized in the case of the detection of an impulsive force, that is, with a time 
duration a much shorter than the mechanical relaxation time (in the absence of feedback), a <C l/7 m . I n fact, 
one could use the following nonstationary strategy: prepare at t = the mirror mode in the stationary state 
cooled by feedback, then suddenly turn off the feedback loop and perform the spectral measurement in the 
presence of the impulsive force for a time T m , such that a <C T m <C l/~f m . In such a way, the force spectrum 
is still well reproduced, and the mechanical susceptibility is the one without feedback (even though modified 
by the short measurement time T m <C l/^ m ). At the same time, the mechanical mode is far from equilibrium 
during the whole measurement, and its noise spectrum is different from the stationary form of Eq. (39), being 
mostly determined by the cooled initial state. As long as T m <C 7m, heating, that is, the approach to the hotter 
equilibrium without feedback, will not affect and increase too much the noise spectrum. Therefore, one expects 
that as long as the measurement time is sufficiently short, the SNR for the detection of the impulsive force (which 
has now to be evaluated using the most general expressions (32) and (34)) can be significantly increased by this 
nonstationary strategy. 

It is instructive to evaluate explicitely the nonstationary noise spectrum of Eq. (34) for the above measurement 
strategy. Simple analytical results are obtained by choosing the following filter function Fr m (t) = 0{t)e~ t / 2Tm 
(0(t) is the Heavyside step function), satisfying J dtFT m (t) 2 = T m . Let us consider the cold damping case first. 
Solving the QLE of the system by taking the equilibrium state in the presence of feedback as initial condition, 
one arrives at the following expression for the detected nonstationary noise spectrum 10 



N Q,non(u) = 7m |Xo(^ - i/2T m )\' 



LU 2 + {l/2T m + lm ) 2 . 2 
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(41) 



where xo is the Fourier transform of the susceptibility in the absence of feedback (see Eq. (29) with g\ = or 
Eq. (30) with g2 = 0) and we have used the high temperature approximation coth(huj/2kBT) ~ 2kBT/tko for 
the Brownian noise. Moreover, (Q 2 ) s t and (P 2 ) s t are the stationary variances in the presence of feedback, 



-oo 2lT 
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(42) 
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Figure 3. Nonstationary noise spectrum for different values of the measurement time, ■y m T m = 10 -1 (dotted line), 
JmTm — 10 -2 (full line), 7 m T m = 10~ 3 (dashed line), j m T m = 10 -4 (dot-dashed line). The figure refers to the cold 
damping feedback scheme, but the curves are indistinguishable from that obtained with the momentum feedback, using 
the same parameters, Q = 10 4 , ( = 10, gi = 32 = 10 3 , k B T/huJm = 10 5 , rj = 0.8. 

where Nqz cd (ui) is given by Eq. (38). The nonstationary noise spectrum for the momentum feedback case is 
analogous to that of Eq. (41), except that one has to use the corresponding stationary values as initial conditions 
(there is an additional term due to the fact that (QP + PQ) s t ^ for momentum feedback 10 ). However, as 
for the stationary case, one can check that the two feedback schemes give indistinguishable results in a large 
parameter region. Therefore we shall discuss the cold damping case only from now on, even though the same 
results also apply to the momentum feedback case with the replacement g 2 ~* gi- 
lt is easy to check from Eq. (41) that the stationary noise spectrum corresponding to the situation with no 
feedback is recovered in the limit of large T m , as expected, when the terms inversely proportional to 7 m T m and 
depending on the initial conditions become negligible, and xq(u> — i / 2T m ) — > Xo(w)- In the opposite limit of small 
T m instead, the terms associated to the cooled, initial conditions are dominant, and since the stationary terms are 
still small, this means having a reduced, nonstationary noise spectrum. This is clearly visible in Fig. 3, where the 
nonstationary noise spectrum is plotted for different values of the measurement time T m , 7 TO T TO = 10 _1 (dotted 
line), 7 m T! m = 10~ 2 (full line), j m T m = 10~ 3 (dashed line), 7 m T m = 10~ 4 (dot-dashed line). The resonance 
peak is significantly suppressed for decreasing T m , even if it is simultaneously widened, so that one can even 
have a slight increase of noise out of resonance. 

The effect of the terms depending upon the feedback-cooled initial conditions on the nonstationary noise 
is shown in Fig. 4, where the noise spectrum is plotted for different values of the feedback gain at a fixed 
value of T m . In Fig. 4a, Nq(uj) is plotted at j m T m = 10~ 3 for 32 = 1 (full line), g 2 — 10 (dotted line), 
g 2 = 10 2 (dashed), g 2 — 10 3 (dot-dashed). For this low value of 7 m T m , the noise terms depending on the initial 
conditions are dominant, and increasing the feedback gain implies reducing the initial variances, and therefore an 
approximately uniform noise suppression at all frequencies. In Fig. 4b, Nq(u>) is instead plotted at j m T m = 1CP 1 
for g 2 = 1 (full line), g 2 = 10 (dotted line), g 2 = 10 2 (dashed), g 2 = 10 3 (dot-dashed). In this case, the feedback- 
gain- independent, stationary terms become important, and the effect of feedback on the noise spectrum becomes 
negligible. 

The significant noise reduction attainable at short measurement times jm,T m -C 1 is not only due to the 
feedback-cooled initial conditions, but it is also caused by the effective reduction of the mechanical susceptibility 
given by the short measurement time, xo{uj) — > xo(<-^ — i/2T m ). This lowered susceptibility yields a simultaneous 
reduction of the signal at small measurement times "f m T m <C 1, and therefore the behavior of the nonstationary 
SNR may be nontrivial. However, one expects that impulsive forces at least can be satisfactorily detected using 
a short measurement time, because the noise can be kept very small and the corresponding sensitivity increased. 
Let us check this fact considering the case of the impulsive force 



}{t) - / exp [-(* - i 1 ) 2 /2 ( x 2 ] cos (uft) , 



(44) 
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Figure 4. Nonstationary noise spectrum for different values of the feedback gain, — 1 (full line), <?2 = 10 (dotted line), 
§2 = 10 2 (dashed), 32 = 10 3 (dot-dashed), with fixed measurement time, 7 m T m = 10 -3 (a), and 7 m T m = 10 _1 (b). (a) 
corresponds to a strongly nonstationary condition, in which the noise is significantly suppressed, thanks to the cooled 
initial condition. In (b) the stationary terms becomes important and the noise reduction due to feedback cooling is less 
significant. The other parameters are Q — 10 4 , ( = 10, hsT/hiVm = 10 5 , n = 0.8. 

where a is the force duration, t\ its "arrival time", and ujf its carrier frequency. The corresponding SNR is 
obtained dividing the signal of Eq. (31) by the nonstationary noise spectra of Eq. (41), and it is shown in Fig. 5. 
As anticipated, the sensitivity of the optomechanical device is improved using feedback in a nonstationary way. 
In Fig. 5a, the spectral SNR, 1Z(ui), is plotted for different values of feedback gain and measurement time. The 
full line refers to g\ = g% = g = 2 • 10 3 and j m T m — 10~ 3 , the dashed line to the situation with no feedback and 
the same measurement time, 5 = and 7 m T m = 10~ 3 ; finally the dotted line refers to a "standard" measurement, 
that is, no feedback and a stationary measurement, with a long measurement time, 7 m T m = 10. The proposed 
nonstationary measurement scheme, "cool and measure" , gives the highest sensitivity. This is confirmed also 
by Fig. 5b, where the SNR at resonance, lZ(oj m ), when feedback cooling is used with g = 2 ■ 10 3 (full line), 
and without feedback cooling (dotted line), is plotted as a function of the rescaled measurement time 7 m T' m . 
The preparation of the mirror in the cooled initial state yields a better sensitivity for any measurement time. 
As expected, the SNR in the presence of feedback approaches that without feedback in the stationary limit 
lmT m 3> 1, when the effect of the initial cooling becomes irrelevant. Fig. 5 refer to a resonant (uif — u> m ) 
impulsive force with j m a = 10~ 4 and y m ti = 3 • 10 -4 , while the other parameters are Q = 10 5 , ^ = 10, 77 = 0.8, 
k B T/huj m = 10 5 . 

The proposed nonstationary strategy can be straightforwardly applied whenever the "arrival time" t± of the 
impulsive force is known: feedback has to be turned off just before the arrival of the force. However, the scheme 
can be easily adapted also to the case of an impulsive force with an unknown arrival time, as for example, that 
of a gravitational wave passing through an interferometer. In this case it is convenient to repeat the process 
many times, i.e., subject the oscillator to cooling- heating cycles. Feedback is turned off for a time T m during 
which the spectral measurement is performed and the oscillator starts heating up. Then feedback is turned on 
and the oscillator is cooled, and then the process is iterated. This cyclic cooling strategy improves the sensitivity 
of gravitational wave detection provided that the cooling time T coo i, which is of the order of 1/ [7m (1 + <?•»)]) 
is much smaller than T m , which is verified at sufficiently large gains. Cyclic cooling has been proposed, in 
a qualitative way, to cool the violin modes of a gravitational waves interferometer in, 6 and its capability of 
improving the high-sensitive detection of impulsive forces has been first shown in. 9 In the case of a random, 
uniformly distributed, arrival time t\ and in the impulsive limit a <C T m , the performance of the cyclic cooling 
scheme is well characterized by a time averaged SNR, i.e., 
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where 1Z(uj, t±) is the nonstationary SNR at a given force arrival time t\ discussed in this section, and TZ(u>, ii) coo z 




Figure 5. (a) Spectrum of the nonstationary SNR, TZ(u), with and without feedback cooling of the initial state. The full 
line refers to a nonstationary measurement, ■y m T m = 10 -3 , in the presence of feedback, g = 2-10 3 (the two feedback schemes 
give indistinguishable curves); the dashed line refers to the no- feedback case, and with the same, short, measurement time 
7m T-m = 10 -3 . Finally, the dotted line refers to a "standard measurement", without feedback, and in the stationary 
limit 7 m T m = 10. (b) Nonstationary SNR at resonance, lZ(u) m ), with and without feedback cooling of the initial state, 
plotted as a function of the rescaled measurement time 7 m T m . The full line refers to the case with feedback-cooled initial 
conditions (g = 2 ■ 10 3 ). The dotted line refers to the no-feedback case, g = 0. The other parameters are uif — u m , 
lm a = 1(T 4 , j m ti = 3 ■ 1(T 4 , Q = 10 5 , < = 10, rj = 0.8, k B T/huj m = 10 5 . 



is the nonstationary SNR one has during the cooling cycle, which means with feedback turned on and with 
uncooled initial conditions. It is easy to understand that 1Z(ui, ti) coo i <C TZ(u>, t%), and, since it is also T coo i <C T m , 
the second term in Eq. 45) can be neglected, so that, 9 

(W> * 7z * dt x n{u>, t x ), (46) 

J-m ~r -L cool Jo 

This time-averaged SNR can be significantly improved by cyclic cooling, as it is shown in Fig. 6, where 
(TZ{lo)) is plotted both with and without feedback. The full line describes the time-averaged SNR subject to 
cyclic feedback-cooling with g = 2 • 10 3 , 7 m T m = 1CP 3 , and T coo i = 1CP 3 T TO . In the absence of feedback, 
in the case of an impulsive force with unknown arrival time and duration er, the best strategy is to perform 
repeated measurements of duration T m without any cooling stage. The measurement time T m can be optimized 
considering that it has to be longer than er, and at the same time it has not to be too long, in order to have a 
good SNR (see the dotted line in Fig. 5)b. In this case, the time-averaged SNR can be written as 

(fcoM) - J" / <ftifto(w,ti), (47) 

where TZa(cu,ti) is the SNR evaluated for g = 0. The dashed line in Fig. 6 refers to this case without feedback, 
and with 7 TO T m = 10~ 3 . The other parameter values are the same as in Fig. 5 and in this case, cyclic cooling 
provides an improvement at resonance by a factor 16 with respect to the case with no feedback. As suggested in 
Ref., 6 one could use nonstationary cyclic feedback to cool the violin modes in gravitational- wave interferometers, 
which have sharp resonances within the detection band. One expects that single gravitational bursts, having a 
duration smaller than the cooling cycle period, could be detected in this way. 




Figure 6. Time averaged spectral SNR with and without cyclic cooling. The full line refers to cyclic cooling with 
7m 7m = 10 -3 , g = 2 • 10 3 , and T coo i = 10~ 3 T m (the two feedback schemes give indistinguishable curves). The dashed line 
refers to the no-feedback case, with the same measurement time jmT m — 10 -3 (see Eq. (47)). The other parameters are 
Lu f = <j m , 7mf = 10 -4 , Q = 10 5 , C = 10, n = 0.8, k B T/huj m = 10 5 . 
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